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Abstract 

We study the effect of a time-delayed feedback within a generic model for a saddle-node bifur- 
cation on a limit cycle. Without delay the only attractor below this global bifurcation is a stable 
node. Delay renders the phase space infinite-dimensional and creates multistability of periodic 
orbits and the fixed point. Homoclinic bifurcations, period-doubling and saddle-node bifurcations 
of limit cycles are found in accordance with Shilnikov's theorems. 
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1. INTRODUCTION 



Time-delayed feedback was originally proposed in the context of chaos control 



Pyra gas. 



1992| as an alternative to the famous OGY method developed earlier by lOtt et all [1990]. 
The idea was to achieve stabilization of unstable periodic orbits (UPOs) by adding, to a 
chaotic system, a control force in the form of the difference between a system variable at time 
t and at a delayed time t — r. This method proved to be very powerful and has been success- 



fully applied to various physica 



was improved by 



2002 



Beck et al. 



Socolar et al. 



200 



2 



Kittel et al. 



syst ems since then Scholl fc Schuster! . 120071 ] . The scheme 



1994 1 , and other variants have been elaborated 



1995 



and applied also t o stochastic systems 



Schlesner et al. 



Goldobin et al. 



2003 



Janson et al. 



20031: 



20041 ] . Moreover, elegant analytical theories 



Baba et al. 



Unkelbach et al 



Hauschildt et al. 



Just et al. 



2003], 



2006 



19971 ] were devel- 



oped supporting, thus, numerical finding s. Apart from the pra ctically relevant application 
of time-delayed feedback, e.g. to lasers [Schikora et all 120061 ] . the interest lies highly on 
the mathematical aspect of the problem. Delay differential equations are difficult to han- 
dle. The delay renders the system infinite-dimensional and the interplay with nonlinearity 
uncovers complex dy namic behaviou r. Delay-induced multistability was already predicted 
in the first paper by iPyragad 19921 ] . The idea that time-delayed feedback may not only 
be used for controlling a system but also for creating new dynamics is not new. Neverthe- 
ess, the investigation of delay-induced bifurcations and multistability is still a growing field 



Balanov et al. 



laser equations 



200 



Xu k Yu 



200 4j with applications to the logistic map as well as to 



Martmez-Zerega &: Pisarchiki . 



2005 : 



Martmez-Zerega et al. 



2003 ] 



2. THE MODEL 



Here we study the in fluence of time-delayed feedback i n a g e neric model representativ e 



for excitable dynamics 



Ditzinger et al. 



1994 



Hu et al. 



1993 



Rappel fc Strogatzj , Il994 | 



The behaviour of the system without delay is governed by a global bifurcation, namely a 
saddle-node bifurcation on a limit cycle {saddle-node infinite period bifurcation, SNIPER), 
that takes place when a certain parameter exceeds a t hreshold. Such a b ifurcation was 



Peinke et al. 



3l 



19891, and also en- 



first observed experimentally in a semiconductor device 
countered in various semiconductor models, e.g. for Gunn domains [Shiau &: Cheng . Il996| . 
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superlattices 


Hizanidis et al.. 


2006; 


Patra et al. 


Krauskopf et al. 


2003; 


Wieczorek et all. 


2002 


1- 



1998] . or lasers 



Dubbeldam et al. 



1999 



20021 ] . Excitability, naturally, yields the system 



highly sensitive to fluctu ations and therefore this s ystem has served perfectly as a n example 



or coherence resonance 



1993] and 



Pikovskv fc Kurths 



Ditzinger et al 



19971 ] shown in the seminal paper by 



Hu et al. 



1994J over ten years ago. In this paper we do not consider the ef- 
fect of noise but extend the generic model by incorporating time-delayed feedback according 
to the Pyragas scheme. The equations are the following: 



x = x(l — x 2 — y 2 ) + y(x — b) — K(x — x r ) (1) 
V = 2/(1 - x 2 - y 2 ) - x(x - b) - K(y - y T ). (2) 

Here x and y denote the variables at time t, while x T and y T the delayed ones at time t — r, 
with r and K being the delay and control strength, respectively. This kind of control is 
called diagonal because the control force may be written in the form of a unity matrix. In 
the absence of delay, i. e. K = 0, b plays the role of the bifurcation paramater. In polar 
coordinates x = r cos tp, y = r sin ip Eq. (JTJ with K = reads 



r = r (l — r 2 ) (3) 
tp = b — r cos ip (4) 

When b < 1 there are three fixed points: an unstable focus at the origin and a pair of 
a saddle-point and a stable node on the unit circle with coordinates (b, +yl — b 2 ) and 
(b, — yl — b 2 ), respectively. The latter collide at b = 1 and a limit cycle r = 1 is born. 
Above but close to the bifurcation, the frequency / of this limit cycle obeys a characteristic 
square-root scaling law / ~ (b — l) 1 / 2 . 



3. LINEAR STABILITY ANALYSIS OF THE DELAY EQUATION 

We prepare the system slightly below the bifurcation (b = 0.95) and switch on the control. 
The first question that arises concerns the stability of the three fixed points and how this 
changes, or not, due to delay. For this, we perform a standard linear stability analysis and 
derive the characteristic equation for the roots, A, which determine the stability of the fixed 
points. For the unstable focus, the characteristic equations is: 
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(1-K + Ke- Ar -Kf + b 2 = 0. (5) 

Due to the presence of the delay, Eq. (J5J) has infinitely many solutions. However, the 
stability of the fixed points is determined by a finite number of critical roots with largest 
real parts. Using the La mbert function W, which is de fined as the inverse function of 
ze z for complex z 



9\?) 
expressed as 



Corless et al. 



19961 : iHald . [I97l|, the solution of Eq. © can be 



A = -W[Kre T{K - 1±b) ] -K + l^ib. 



(6) 



In the case of the saddle and the node, the characteristic equation can be factorized into 
two equations: 



A + K + 2- Ke- Ar = (7) 
A s ' n + K =f - Ke~ AS ' nr = 0, (8) 



with solutions: 



Ai = -W[Kre ri2+K) } — 2 — K (9) 

T 

A2 s,n = \ W [ KTe <KTVT^)^ _ K± ^1^1?. (10) 
T 

The superscripts 's' and V denote the saddle (upper sign) and the node (lower sign), 
respectively. Figure [1] shows the real parts of the eigenvalues A as a function of r for a 
fixed value of K for all three fixed points. One may see the eigenvalues of the uncontrolled 
system at r = 0, and their interaction with the delay- induced modes (blue) with increasing 
t. In all three cases, control is unable to change the stability of the fixed point: in the 
case of the unstable focus, the mode with the largest real part (red) tends to very small 
values with increasing delay, remaining however positive. The same holds for the unstable 
mode of the saddle. Symmetric behaviour is observed for the stable modes of both saddle 
and node: they tend to zero as a function of r but remain negative. The picture does not 
change qualitatively even for other values of K and therefore one might conclude that no 
delay-induced bifurcations of fixed points take place. 
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FIG. 1. Real parts of the complex eigenvalues A as a function of r, for fixed K = 1 and b — 0.95. For the 
(a) unstable focus, (b) the saddle point and (c) the stable node. The modes emerging from the 
uncontrolled system and the delay-induced modes are marked red and blue, respectively. 

4. GLOBAL BIFURCATION ANALYSIS 

However, the above local analysis gives no information on the global changes in phase 
space that delay potentially induces. A numerical investigation shows, in fact, that there 
exists bistability in a certain parameter regime in the K-r plane: trajectories starting close 
to the saddle point are attracted by a delay-induced limit cycle, whereas trajectories starting 
elsewhere end up in the stable node. Keeping r = 3 fixed we find the critical value K c of K 
for which this delay-induced limit cycle is born and observe a scaling T ~ In \ K — K c \ in the 
period T of the corresponding oscillations, typical for the case of a homoclinic bifurcation. 
In Fig. [2] phase portraits of the system below, at and above the bifurcation are shown. 
Trajectories with different initial conditions are shown: one starting from the vicinity of the 
unstable focus in the origin (blue) and one from the vicinity of the saddle (red). Bistability 
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is revealed in Fig. [2^c) where two attractors (the stable node and the delay-induced limit 
cycle) coexist. The 'kink' in the trajectory shortly before the loop closes is due to the control: 
the control force starts acting at t — 3 when the system is still moving on the slow part of 
the unit circle. Therefore, its effect is not so noticeable. As the system moves faster the 
control force attains higher values and the trajectory starts deviating from its deterministic 
path at t = 13. This deviation becomes large at t = 18 where the trajectory, shortly before 
settling in the stable node, appears to be "attracted" to the saddle, resulting in this 'kink' 
in the x — y projection. Also, as K approaches the critical value, the trajectory passes closer 
and closer to the saddle on its way to the stable node. This ends in a homoclinic orbit at 
K = K c , (Fig. EJjb)) from which a periodic orbit is generated (Fig. EJ^c)). The period of the 
born limit cycle scales according to T ~ — A u ~ 1 ln\K — K c \, where A n is the real part of the 
least unstable eigenvalue of the saddle point (i. e. , the one closest to the imaginary axis). 
One may calculate this from Eq. (ED) and find A u _1 = 0.1739 -1 = 5.75. This is in rather 
good agreement with the slope of the solid line in Fig. [2](d) which equals 5.35. 



In the following we use a bifurcation continuation tool Engelborghs et all 120021 . 12001 ] 



and follow the homoclinic bifurcation in the K — r plane. The produced bifurcation curve 
can be seen in Fig. [3] (left). It consists of two main curves: one running through points A — E 
and a second tongue-like curve. In the white area the system is monostable (stable fixed 
point) while in the yellow area a delay-induced periodic attractor is born via a homoclinc 
bifurcation marked by the red curves. 

At this point one should emphasize the role of the saddle-point: due to the delay, the 
saddle possesses no longer two distinct eigenvalues (one positive, i. e. unstable, and one 
negative, i. e. stable) but infinitely many. Moereover, complex eigenvalues come into play 
as well. The eigenvalues, however, that determine the behaviour of the colliding homoclinic 
orbit, are the leading ones, i. e. those closest to the imaginary axis. In Fig. H] one can see 
the eigenvalue spectrum for two parameter values on the homoclinic curve and notice that 
the leading eigenvalues of the saddle (red) are a positive real eigenvalue, as in the original 
uncontrolled system, and a complex conjugate pair with negative real parts, generated by the 
delay. This means that the saddle may turn into a saddle-focus for certain values of K and 
r. Homoclinic orbits attached to a saddle-focus approach the fixed point in an oscillating 
manner. This explains the phase portraits in Fig. [3] (right) which become more and more 
complicated as K increases. 
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FIG. 2. (a) Two dimensional projection of the phase space below the homoclinic bifurcation (K = 0.335). 
(b) Homoclinic orbit (red) achieved at K c — 0.3401. (c) Delay-induced limit cycle (red) above the 
homoclinic bifurcation (K = 0.3438). (d) Scaling of the oscillation period T above but close to the critical 
point K c (crosses: simulation data, solid line: linear fit). Full and open circles mark stable and unstable 
fixed points, respectively. Parameters: b — 0.95, r = 3. 



From the above, it is clear that the two basic ingredients responsible for the delay- 
induced dynamics in our system are the homoclinic orbits and the saddle-foci. The the- 



ory of homoclinic bifurca t ions 



Guckenheimer fc Holmesl . 



1986 



or ordinary differential equations has been well developed 



Kuznetsov 



1995 



Wiggins! . 



19881 ] . As already mentioned, 
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FIG. 3. Curve of homoclinic bifurcations (red) in the K — r plane (left). A - E labels various points with 
homoclinic orbits, which are shown in the x — y phase plane in the panel on the right. Delay- induced limit 
cycles exist, in addition to the stable fixed point, in the yellow area. The blue dashed curve separates the 
regions ctq < (left) and ctq > (right) 



global bifurcations are strongly related to excitability and therefore one expects to encounter 
them in excitable syste ms. Various physical syst ems such as modulatio n-doped semicon- 



ductor heterostructures 



2003 



Wieczorek et al 



Bordyugov fe Engel 



200jJ have been studied in this respect, both theoretically and ex- 



20021, neuron models 



Dottling fe Scholll . Il992fl. semicond u ctor lasers [Krauskopf et al. 



Feudel et al. 



2000] and chemical systems 



perimentally. On the other hand, less work has been carr ied out for systems with delay 



undergoing such nonlocal bifurcations 



Sethia fe Sen! . 120061 ] . It is therefore appropriate to 



analyze a generic system like the one studied here which, despite its simplicity, exhibits rich 
delay- induced dynamics with a homoclinic bifurcation as key component. 
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FIG. 4. Spectra of the eigenvalues A of the saddle for two points on the homoclinic bifurcation line: (a) 
K = 0.3401, r=3 (B in Fig. EJ and (b) K = 0.57, r = 7.28 (E in Fig. [3]). The saddle quantity 
o o = Re(A s ) + Re(A u ) is negative in (a) whereas in (b) it is positive, (b = 0.95) 



DELAY-INDUCED MULTISTABILITY 



In what follows we will app ly the theorems o n homoclinic orbits connecting to saddle 



foci as developed by Shilnikov 



Kuznetsov 



1995]. According to them, the so-salled saddle 



quantity is crucial for the homoclinic bifurcations occurring in high-dimensional systems. 
The saddle quantity is defined as o"o = Re(A s ) + Re(A u ), where A s and A u are the leading 
stable and unstable eigenvalues, respectively. Shilnikov proved, among other, that negative 
cr results in the birth of a unique stable limit cycle from a homoclinic orbit. On the other 
hand, for ctq > 0, a wide variety of homoclinic bifurcations may occur, some of which involve 
infinitely many periodic orbits in the vicinity of the homoclinic orbit. The green dashed 
curve in Fig. [3J shows the condition = 0. Along the homoclinic bifurcation line in Fig. [3J 
the saddle quantity changes sign, thereby allowing for both scenarios to take place. Figure 0] 
shows the eigenvalue spectra for two different points on the bifurcation line corresponding 
to negative and positive saddle quantities, respectively. 

In the following, we restrict ourselves to a fixed value of r = 7 and reveal multistability 
beyond the homoclinic bifurcation. For r = 7, a homoclinic bifurcation takes place at 
K = 0.17145 (point A in Fig. [3D- In this case cro = —0.0116, and the bifurcation creates 
one stable limit cycle, with Floquet multipliers within the unit circle (inset of Fig. [5]^a) 
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FIG. 5. (a) Period T of limit cycle born in a homoclinic bifurcation at (if, r) = (0.17145, 7) (point A in 
Fig. O (T < 0). (b) Period T of limit cycles in the multistable regime at (if, r) = (0.213, 7) (point Z in 
Fig. O Co > 0), undergoing infinitely many fold (F) and period-doubling (PD) bifurcations, before ending 
in a homoclinic orbit h for T — > oo at if = 0.213. Solid blue and red dashed lines denote stable and 
unstable limit cycles, respectively. The insets show the two leading Floquet multipliers of the periodic 
orbit /ii — 1 (green) and (12 (blue) with if as a parameter in (a), and T as a parameter in (b), in the 
complex plane, b = 0.95 



). The period T of this limit cycle increases monotonically as the bifurcation point is 
approached (Fig. G^a)). Moving further along r = 7 other homoclinic bifurcation curves are 
crossed, e.g. at K = 0.213 (Fig. [5(b), cf. point Z in Fig. [3]). There, the saddle quantity is 
positive (ctq = 0.0023, calculated analytically from Eq. (Till ), and the picture is much more 
complicated: an infinite number of bifurcations take place, which are related to saddle-node 
(fold) bifurcations of pairs of stable and unstable limit cycles, and additional period-doubling 
(flip) bifurcations of the stable limit cycles. The dependence of the period T of the limit 
cycles upon if is a nonmonotonic multivalued function, whose turning points are associated 
with saddle-node bifurcations. In between the fold bifurcations of the stable limit cycles, 
pairs of forward and inverse period doubling bifurcations occur. In the insets of Fig. [5] 
the trivial Floquet multiplier //} = 1 and the leading Floquet multiplier /x 2 of the periodic 
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orbit are plotted. It can be seen how [12 changes as T is varied along the multivalued 
function in the main figure, showing how fold and flip bifurcations occur, at /i± = 1 and 
/ii = —1, respectively. One should also expect other bifurcations near the critical point due 
to secondary homoclinic orbits which are beyond the scope of this paper. 



6. CONCLUSIONS 



In conclusion, we have presented a mechanism for delay-induced multistability in a system 
near a global bifurcation. In addition to the fixed point attractor which the uncontrolled 
system already possesses, a time-delayed feedback in the form of Pyragas difference control 
induces one or more coexisting limit cycle attractors. Depending upon the feedback control 
strength K and the delay time r, either a single stable limit cycle is born in a homoclinic 
global bifurcation, or an infinite number of (stable and unstable) periodic orbits is induced 
undergoing a rich menagerie of bifurcation scenarios including period doubling and fold 
bifurcations. We have shown that the key ingredient in the observed dynamics is a homoclinic 
orbit connected to a saddle-focus created by delay. A bifurcation continuation in the K — t 
plane was performed. Moreover, we were able to verify Shilnikov's theory of homoclinic 
bifurcations in a certain parameter regime. The excitable nature of the system and the 
infinite-dimensional phase space, due to delay, appear to play a crucial role in the induced 
homoclinicity. 

These results are interesting also from the point of view of applications, since our generic 
model is representative for a wide range of real-world systems. For instance, the transition 
from stationary to moving field domains in semiconductor superlattices has been shown 
to be a ssociated with a sadd le-node bifurcation on a limit cycle as described by Eq.flU) at 



K = 
system 



iizanidis et al 



Schlesner et al. 



2006], and time-delayed feedback control can also be realized in this 



20031 ] . Already without delay, this system has been noted for its 



high multistability of stationary domain states 



Kastrup et al. 



199 



1996 



Prengel et al. 



1994j, and bistability or 



nanostructures, see e.g. 



higher multistability 



Meixner et al. 



ras be e n fou nd in many other semiconductor 



2000j: 



Scholll [2001 1. 
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